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ABSTRACT 



We first consider the Lagrangian formulation of general relativity for perturbations 
with respect to a background spacetime. We show that by combining Ncether's method 



> 

in 

with Belinfante's "symmetrization" procedure we obtain conserved vectors that are inde- 
pendent of any divergence added to the perturbed Hilbert Lagrangian. We also show that 
the corresponding perturbed energy-momentum tensor is symmetrical and divergenceless 



On 

but only on backgrounds that are "Einstein spaces" in the sense of A.Z. Petrov. de Sitter 
or anti-de Sitter and Einstein "spacetimes" are Einstein spaces but in general Friedmann- 
Robertson- Walker spacetimes are not. Each conserved vector is a divergence of an anti- 
symmetric tensor, a "superpotential" . We find superpotentials which are a generalization 
of Papapetrou's superpotential and are rigorously linear, even for large perturbations, in 
terms of the inverse metric density components and their first order derivatives. The super- 
potentials give correct globally conserved quantities at spatial infinity. They resemble Ab- 
bott and Deser's superpotential, but give correctly the Bondi-Sachs total four-momentum 
at null infinity. 



* A summary of this work has been presented at the Conference on "Fundamental 

Interactions: From Symmetries to Black Holes" , Brussels, March 25-27, 1999. It has been 

published in the Proceedings of the Conference under the title Conservation Laws for 

Large Perturbations on Curved Backgrounds with ref.: Fundamental Interactions: From 

Symmetries to Black Holes (eds.: J.M. J.M. Frere, M. Henneaux, A.Servin & Ph. Spindel), 

p.p. 147-157, Brussels: Universite Libre de Bruxelles (1999), see also ^r-qc/9905088 . 
f E-mail: petrov@sai.msu.su 
X E-mail: jkatz@vms.huji.ac.il 



1 



Next we calculate conserved vectors and superpotentials for perturbations of a Fried- 
mann- Robertson- Walker background associated with its 15 conformal Killing vectors given 
in a convenient form. The integral of each conserved vector in a finite volume V at a given 
conformal time is equal to a surface integral on the boundary of V of the superpotential. 
For given boundary conditions each such integral is part of a flux whose total through a 
closed hypersurface is equal to zero. For given boundary conditions on V, the integral 
can be considered as an "integral constraint" on data in the volume and this data always 
includes the energy-momentum perturbations. We give explicitly these 15 integral con- 
straints and add some simple applications of interest in cosmology. Of particular interest 
are Traschen integral constraints in which the volume integral contains only the matter 
energy-momentum tensor perturbations and not the field perturbations. We show that 
these particular integral constraints are associate with time dependent linear combina- 
tions of conformal Killing vectors. Such linear combinations are neither Killing vectors nor 
conformal Killing vectors. We also find that if we add the "uniform Hubble constant hy- 
persurface" gauge condition of Bardeen, there exists 14 such integral constraints. The 
exception is associated with conformal time translations (k = ±1) or conformal time ac- 
celerations (k = 0). As an example we find the constants of motion of a spacetime that is 
asymptotically Schwarzschild-de Sitter (k = 0). 
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1. Introduction 

(i) Conservation laws and cosmology. 

Conservation laws associated with "symmetric" infinitesimal displacements in Fried- 
man-Robertson-Walker spacetimes have been used in relativistic cosmology on several 
occasions. Infinitesimal displacements are characterized by vector fields and, as we shall 
see, the vectors used in some applications are not always Killing vectors nor even conformal 
Killing vectors. 

An example in which no Killing nor conformal Killing vectors are used has been given 
by Traschen [I] who introduced "integral constraints" in terms of "integral constraints 
vectors" . Traschen and Eardley [2] analyzed measurable effects of the cosmic background 
radiation due to spatially localized perturbations. By using "integral constraints" they 
pointed to an important reduction of the Sachs- Wolfe [3] effect on the mean square angular 
fluctuations at large angles of the cosmic background temperature due to local inhomo- 
geneities. Traschen's integral constraints vectors have a somewhat intriguing origin [4]. 
The equations for integral constraint vectors have been studied by Tod [5]. He showed 
that these equations are conditions for a spacelike hypersurface to be embeddable in a 
spacetime with constant curvature of which the solutions are Killing vectors. In Katz, 
Bicak and Lynden-Bell [6], a paper referred to as KBL97, integral constraints appear as 
conservation laws with Killing vectors in a de Sitter background; more on this below. 

Local differential conservation laws, rather than global ones, have been used by Veer- 
araghavan and Stebbin [7] . They found and used a conserved "energy- momentum" pseudo- 
tensor in an effort to integrate Einstein's equations with scalar perturbations and topo- 
logical defects in the limit of long wavelengths on a Friedmann-Robertson- Walker space- 
time with flat spatial sections (t = const, k = 0). Uzan, Deruelle and Turok [8] realized 
that these conservation laws might be associated with the conformal Killing vector of 
time translations and they extended Veeraraghavan and Stebbin's method to Friedmann- 
Robertson- Walker perturbed spactimes with non-flat spatial sections (k = ±1). More on 
this in section 4. 

In Lynden-Bell, Katz and Bicak's [9] study of Mach's principle from the relativistic 
constraint equations, conservation laws yield a general proof that the total angular mo- 
mentum (and the total of any conserved perturbation of the current which deriving from 
a "superpotential" ) must be zero in any closed universe. 

As a final example we mention KBL97's analysis of the globally conserved quantities 
that result from mapping a Friedmann-Robertson- Walker perturbed spacetime on a de 
Sitter background with its ten Killing vectors. 

With these different examples in mind, it made good sense to study the properties 
and physical interpretation of conservation laws and their superpotentials, in the context 
of relativistic cosmology, associated with arbitrary displacements in a background as was 
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done in KBL97. In fact, the theory has wider applicability than relativistic cosmology 
since the background may be any spacetime and there are plenty of examples in general 
relativity in which backgrounds are used, 
(ii) Noether's method and its problems. 

KBL97 used the fairly standard method of Ncether (see for instance Landau and 
Lifshitz [10]) to derive conservation laws from the Lagrangian Cq of the perturbations of 
the gravitational field** 

c G = c-I, c = -±-(R + d^), I = -±-(R + d^), « = ^- (i.i) 

Here a * means multiplication by y/—g, a bar referees to the background, R is the scalar 
curvature and /c M is some vector density. Noether's method associates a conserved vector 
density I M with any vector £^ that generates small displacements. It applies to pertur- 
bation theory on any background and provides a "canonical" energy-momentum tensor 
perturbation. Moreover, the conserved vector is always [11] the divergence of an antisym- 
metric tensor, a superpotential, (1^ = d v I^ v ', = —I u ^). This has the great practical 
advantage that integrals of complicated 7 M 's in a volume V are equal to often much sim- 
pler integrals of J Mi "s on the boundary of V. Noether's method is the most direct and 
easy way to construct superpotentials, field energy tensors and conserved vector densities 
with arbitrary backgrounds and for arbitrary £^ though the same results can of course be 
worked out directly from the perturbed Einstein equations. But, at least in our case, this 
is far more complicated than with the method developed here as we shall see. 

Noether's method has, however, two unsatisfactory features. First the Lagrangian 
density is not unique. A divergence <9 M /c M can and must be added to the Hilbert Lagrangian 
because the latter leads [12] to Komar's [13] conservation law which gives the wrong mass 
to angular momentum ratio with an "anomalous" factor of two in the weak field limit [14] 
and does not give the Bondi mass [15] at null infinity [16]. Divergences are also added to 
comply with different boundary conditions. Various divergences have thus been added to 
R for different reasons. M0ller [17], using a tetrad representation with g^ v = rjape^e^, 
would have taken a /c M = g^ v ' g po r (e"D ff e^-e"D pC^Va/i- York [18], using a foliation, would 
have chosen = 2 (en^D u h u — n u D u h^) with (n^n^ = e = ±1) the normal vectors 
of his closed hypersurfaces. KBL97 wanted a field energy tensor quadratic in first order 
derivatives and took therefore like Rosen [19] a long time before 

& = -^=D v {-ggr). (1-2) 



** Notations are properly defined in section 2. Here we assume the reader to be familiar 
with current notations. 
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_D M is a covariant derivation with respect to the background metric g^ v . Second, the 
canonical field energy momentum tensor is not symmetrical nor is it divergenceless. On a 
flat background, the energy-momentum tensor is divergenceless but is still not symmetrical 
and the angular momentum is not conserved; it does not include the helicity of the field. It 
thus appears that conservation laws obtained by Ncether's method have an unsatisfactory 
weak field limit on a flat background at least as far as angular momentum is concerned. 

To remedy that situation we suggest in this paper to modify Ncether conserved vectors 
using Belinfante's [20] trick in classical field theory. It is an easy matter to adapt his 
method to perturbation theory on curved backgrounds. Belinfante's modification leads to 
energy-momentum tensors which ensures, at least in classical field theory, that angular 
momentum includes the helicity and is then conserved. The Belinfante trick has been 
applied by Papapetrou [21] to general relativity in an effort to calculate the total angular 
momentum at spatial infinity. 

These new conserved vectors have none of the drawbacks just described and in addition 
have very appealing new properties, 
(iii) A summary of theoretical results. 

It may be useful, at this stage, to give a summary of our main theoretical results, we 
mean those valid in general, not only in relativistic cosmology. 

(a) We find that there exist a conserved vector density X M associated with any vector 
generating infinitesimal displacements and which is the divergence of a superpotential 
T^ v ; it is of the following form 

= r*i v + Z» = d v T? v . (1.3) 

q-ji re p resen ts a matter plus field energy-momentum tensor density perturbation relative 
to the background. is a vector density that is only equal to zero if is a Killing vector 
of the background, which we denote then by As a consequence of Eq. (1.3) a volume 
integral F of X^ equal the surface integral of Z^ v on its boundary S 

F= f Z^dV^ = I l^dS^ (1.4) 
Jv Js 

and the total flux F through a closed hypersurface V is equal to zero. That is what is 
meant X^ being a conserved current. Equalities such as (1.4) may be regarded as "integral 
constraints" in the following sense. Suppose that boundary values on S and thus Z^ v 
are known. Then Eq. (1.4) represents constraints on the perturbations of the energy- 
momentum tensor which is always part of TJ\ There are thus as many integral constraints 
as there are displacement vectors £^ (but not all are equally interesting). Our notion 
of "integral constraints" is slightly different from that introduced by Traschen [1]. She 
calls integral constraint an expression like Eq. (1.10) below in which the boundary term 
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is equal to zero and the volume integral contains only the perturbation of the matter 
energy-momentum tensor. One may wish to see in Eq. (1.4) a definition of quasi- local 
conservation laws if the boundary is not at the border of spacetime itself. 

(b) The conserved vector, and thus the corresponding are independent of any divergence 
added to the Hilbert Lagrangian density of the perturbations. It has been noticed before 
by Bak, Cangemi and Jackiw [22] that Belinfante's modification of the Ncether currents 
obtained from Hilbert's or Einstein's Lagrangians lead to the same symmetric and diver- 
genceless energy-momentum tensor relative to a flat background in Minkowski coordinates. 
Ours is a generalization of this finding for any divergence added to the Hibert Lagrangian, 
for arbitrary perturbations with respect to any background in arbitrary coordinates. 

We want to stress that since our conserved vectors are independent of an added di- 
vergence, they are also independent of boundary conditions. This result is in line with 
classical field ideas. The opposite view that pseudotensors and superpotentials must de- 
pend on boundary conditions has been held for instance in [23]. 

(c) From Eq. (1.3) follows that for each Killing vector of the background there exists a 
conserved vector and a corresponding superpotential such that 

This expression looks also very much like a conservation law in classical field theory. 

(d) = T^g pi/ is symmetrical and divergenceless if and only if the background is an 
Einstein space in the sense of A.Z. Petrov [24], that is if R^ v = —Kg^ where A is necessarily 
a constant, de Sitter spacetimes belong to that category. Other Friedmann-Robertson- 
Walker spacetimes that are currently used in cosmology do not. 

(e) The new superpotential is reminiscent of many well known ones (see section 2 for 
details) and has a simple form: 

= = ip\MD £*] + _D a ( iPlMg*]* - r^g^p] £ (1.6) 

K K V / 

In this expression l pv is the perturbed inverse metric density: 

i^=g^-g~^. (1.7) 

The superpotential has the remarkable property of being linear in . Linearity is a 
valuable property; the linear approximation is not different from the non linear one. Global 
exact conservation quantities of know asymptotic fields with unknown sources can be 
calculated and given physical meaning. The superpotential (1.6) satisfies standard criteria 
of global conservation laws in asymptotically flat spacetimes at spatial and at null infinity 



6 



(see appendix). Also second order corrections of the energy-momentum tensor due to field 
energy contributions are readily calculable from our formulas. 

Those are the principal theoretical results of the paper, 
(iv) The 15 Conformal Killing vectors of cosmological backgrounds and their associated 
conservation laws and integral constraints. 

To illustrate our new conservation laws in theoretical cosmology we consider the con- 
served vectors and superpotentials associated with the 15 conformal Killing vectors of 
Friedmann-Robertson- Walker spacetimes. By conformal Killing vectors we mean the 15 
linearly independent solutions of the conformal Killing vector equations 



The 15 conformal Killing vectors include by definition the 6 pure Killing vectors for which 
D p ^ p = 0. Friedmann-Robertson- Walker spacetimes are conformal to Minkowski's space- 
time. There are thus similarities between the 4 Killing vectors of translations, the 3 rota- 
tions, 3 Lorentz boosts, 3 center of mass position, 1 dilatation and the 4 "accelerations" 
of Minkowski's spacetime. Such similarities are helpful in geometrical interpretations. A 
presentation of those conformal Killing vectors of Minkowski spacetime in a form that ap- 
peals to physicists is given in Fulton, Rohrlich and Witten [25] who have a particular liking 
for the accelerations which they explain well and to which we referee the reader interested 
in those slightly unfamiliar conformal Killing vectors . 

Now follows a brief summary of the results, 
(a) We give the 15 linearly independent solutions of Eq. (1.8) in a simple mathematical 
form. We are not interested in the algebra of the conformal group. We are mostly interested 
here in quasi-local or in global conservation laws or integral constraints for volumes in 
a sphere (parametrized by r) at a given instant of conformal time r\. Thus, in (77, x k ) 
coordinates we are looking for integrals of the form 



Any linear combination of such integrals with time dependent coefficients are still of the 
same form. We give what seem to us the 15 simplest linear combinations. Notice that 
linear combinations of conformal Killing vectors with time dependent coefficients are not 
conformal Killing vectors anymore. Notice also that the left hand integrand of such linear 
combinations is not a zero component of a conserved vector anymore, 
(b) We next turn our attention to those linear combinations in which the volume integrand 
depends only on the matter energy-momemtum perturbations 5T® thus of the form 



(1.8) 




(1.9) 




(1.10) 
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There are 10 integral constraints of this form, 6 are associated with the pure Killing 
vectors of the Friedmann-Robertson- Walker spacetimes and 4 are the vectors found by 
Traschen [1]. Thus Traschen's "integral constraint vectors" appear here as linear combi- 
nations of conformal Killing vectors with time dependent coefficients. 

(c) We then show that if we apply the uniform Hubble expansion gauge studied by Bardeen 
[26] all but one of the 15 conformal Killing vectors are associated with integral conservation 
laws of the form (1.10). The exception is associated with conformal time translations if 
k = ±1 or conformal time accelerations if k = 0. A look at Eq. (1.10) shows that these 
integrals might be constructed directly from Einstein's constraint equations. However the 
V M 's though simple are not all that easy to guess as we shall see. The integral constraints 
have often simple geometrical interpretations by analogy with classical mechanics. With 
conformal Killing vectors they are momenta of order 0, 1 or 2. 

(d) We also give a non-trivial example of globally conserved quantities on a background that 
is not asymptotically flat. We calculate the constants of motion of spacetimes that are 
asymptotically Schwarzschild-de Sitter (with k=0). Because of the high degree of symmetry 
of the asymptotic conditions we find 13 globally conserved quantities equal to zero, that is 
13 Traschen's like integral constraints. There are 2 constants of motion which are different 
from zero. 

Such are the examples of cosmological interest studied in this paper, 
(v) Presentation of the paper. 

In the following section we describe in detail the way to obtain Ncether's conserva- 
tion laws on curved backgrounds. A few parts of that chapter are taken from KBL97 
but they are well worth repeating here. Section 3 introduces the Belinfante correction to 
Ncether's conserved vectors. There we show that the modified conserved vectors and their 
associated superpotentials are unchanged if we add a divergence to the Hilbert Lagrangian 
for the perturbations. At the end of section 3 we derive Rosenfeld's [27] identities which 
give beautiful relations among complicated quantities of interest. Rosenfeld's identities are 
refereed to as "cascade equations" by Julia and Silva [28]. In that section we also obtain 
the theoretical results summarized above in (iii). In section 4 we briefly describe how 
we found the 15 conformal Killing vectors of Friedmann-Robertson- Walker spacetimes in 
appropriate coordinates and how to calculate the corresponding conserved vectors and 
superpotentials in a 1 + 3 standard decomposition. The same section contains the exam- 
ples just described in (iv). In section 5 we emphazise the role of superpotentials in the 
development of conservation law theory for general relativity and make a short review of 
historically important superpotentials on a flat background. A flat background is quite 
useful to study conservation laws in general relativity as pointed out by Rosen [19] a long 
time ago. 
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Finally in appendix we briefly show that our new superpotential has the normally 
expected global properties at null and spatial infinity in asymptotically flat spacetimes. 

Each section is preceded by a summary which gives the motivations and point out 
where the reader will find the principal formulas. The body of the sections themselves 
are written for readers who are interested in the mathematical details. However, most 
elaborate but straightforward calculations are not given in detail. Unfortunately general 
relativity is replete with them. 



2. Noether 

(i) Motivations and summary of results. 

In this section we apply Ncether's method to the Lagrangian Cq defined in Eq. (1.1) 
with the vector density /c M defined in Eq. (1.2). This is the KBL97 Lagrangian. The 
conserved vector I p obtained in this way is given in Eq. (2.17). A look at Eq. (2.17) 
shows that it contains a canonical energy-momentum tensor density 9^ explicitly written 
in Eq. (2.18) with i% given in Eq. (2.13), an "helicity" term a ppcr given by Eq. (2.14) and 
a vector density see Eq. (2.19), defined in terms of the derivatives of the Lie derivatives 
of the background metrics g pa or z p(J defined in Eq. (2.11); rf is thus zero if £^ is a Killing 
vector ^ of the background. 

The canonical energy-momentum tensor density is neither symmetrical nor diver- 
genceless except on a flat background [see Eq. (5.1a)] in which case the canonical field 
energy-momentum i% reduces to Einstein's pseudo-tensor in Minkowski coordinates X^. 
In the same coordinates the helicity tensor density a^ pa ^ is that which was given by Pa- 
papetrou [21]. The conserved vector density I p is thus a generalization of standard old 
results to finite perturbations of a curved background with arbitrary vectors not only with 
Killing vectors. I p is also written as the divergence of a superpotential which is given 
in Eq. (2.21) and also by Eq. (2.22). 

The present section also contains the energy-momentum tensor, the helicity, and the 
superpotential that are obtained from the Hilbert Lagrangian density 

t' G = ~{R-k). (2.1) 

The quantities, 9' p and a ,p,pa are defined by Eqs. (2.24); the superpotential K pv is given 
by "Eq. (2.21) minus its £fc-term". While 0% and, to a certain extend, also <7 MP<J reduce to 
familiar quantities on a flat background, 9 lp and a'^ p ^ are much more complicated; 9'Jf 
contains second order derivatives. In fact we gain in clarity and simplicity by starting the 
calculations with Cq rather than C' G ; the end product in section 3 is independent of <9 M /c M . 

The conservation law d^I p = has been studied in KBL97. Therefore the present 
section is mainly mathematical; it gives the necessary ingredients for sections 3 and 4. 
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(ii) The conserved vector density t M . 

Let g pu (x x ) be the metric of the perturbed spacetime M. and g pu be the metric of the 
background M. both with signature —2. Once we have chosen a smooth global mapping 
such that each point P of Ai is mapped on a point P of Ai, we can use the convention that 
P and P shall always be given the same coordinates — x^ . This convention implies that 
coordinate transformations on M. inevitably induce the same coordinate transformations 
with the same functions on M.. With this convention, such expressions as g^ — g pl , become 
true tensors. However if the particular mapping has been left unspecified, we are still free 
to change it. The form of the equations for perturbations must inevitably contain a gauge 
invariance corresponding to this freedom. 

Le 
follows 



Let R\ pa and R vpa be the curvature tensors of M and M. These are related as 



R\ pa = D p Al - D a At p + A^AL - A^A*, + r\ P(J . (2.2) 

Here D p are covariant derivatives with respect to g pu and A^ is the difference between 
Christoffel symbols in M. and M.: 

= r*, - = \g Xp (D p9pu + D v g pp - D p g pv ) . (2.3) 

Our quadratic Lagrangian density Cq for the gravitational field is here defined by Eq. 
(1.1) with k^ given in Eq. (1.2). The caret means, as we said before, multiplication by 
\/—g, never by \J —~g. Thus, if R = y/—gR, R will unambiguously mean y/ 1 —~g~R. Notice 
that R = yJ—gR ^ R = R. The vector density k^ can also be written in the following 
form that is often useful in calculations: 

& = -^=D V {-g 9 n = g^A* - g p °A» (2.4) 



d p k^ cancels second order derivatives of g pv in R. C is the Lagrangian density used by 
Rosen [19]. C is C in which g pu has been replaced by g pu . When g pv = # , Cq is thus 
identically zero. The following formula, deduced from Eqs. (2.2) and (2.4), shows explicitly 
how Lq is quadratic in the first order derivatives of g pv or, equivalently, quadratic in A^ pa : 

to = tJ" (a - a - " A - A - } " h z> ^- (2 - 5) 

where is the perturbed metric density defined in Eq. (1.7). If the background 
is flat and denoted Mo and if we use Minkowski coordinates then g piy = n pu = 
diag(l, -1,-1,-1), T* v = and 

£g = ~ (R + = (r^ - rp^) (2.6) 
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which is Einstein's [29] Lagrangian. Lq is thus a generalization of Einstein's Lagrangian 
density to perturbations on a curved background. 

Lie differentials are particularly convenient in describing infinitesimal displacements 
in both Ai and Ai; if the mapping was defined before the displacements it remains defined 
after displacements. Let Aa? M = £ M AA represent an infinitisimal one-parameter displace- 
ment generated by a sufficiently smooth vector field £ M , the corresponding infinitesimal 
change in tensors are given in terms of Lie derivatives with respect to this vector field 
^■9fiu =££9lj,isA\, etc. The Lie derivatives may be written in terms of partial derivatives 
<9 M , covariant derivative with respect to g^, or covariant derivative _D M with respect 
to g^ u . Thus, 

£{-9vlv = 9^\d^ X + 9u\d^ X + t X d x g^ u (2.7a) 
= g^xD^ + g„ x D^ x + £ x D x g^ (2.76) 
= g^ x D^ x + g uX D^\ (2.7c) 

Consider now the Lie derivative £^C of t in Eq. (1.1), not of Cq- The Lie derivative of a 
scalar density like t is the ordinary divergence <9 M (££ M ). With the variational principle in 
mind we can thus write the following identity 

£^C = j-G^£^g^ - (g pa £^ pa - g^£^ pa + £^) = d p (££ M ) (2.8) 

where Einstein's tensor density G^ v = R^ v — ^g^R is the variational derivative of 2kC 
with respect to g pv . Equation (2.8) is easily converted into a conservation law by treating 
the first term after the first equality sign as follows: (a) replace £^g [liV by the expression 
(2.7c), (b) use the contracted Bianchi identities D V G^ V = and (c) use Einstein's equations 
= kT£. Identity (2.8) becomes then a conservation law of this form 

= 0, with <" = tft v - ^ (g pa £^% ~ + ~ ( 2 - 9 ) 

Now comes another exercise which consists in replacing £^g fIL , and the D-derivatives of 
£i9^v that appear in £^T^ a - see Eq. (2.20) below - by D-derivatives and D D-derivatives 
using Eq. (2.76) this time. The relation between the two kinds of derivatives, D and D is 
best illustrated on the following simple case 

B v e = d„e + K P i P = + A^p. (2.10) 

The relations is: (a) write D-derivatives in terms of d-derivatives and T's (b) replace the 
<9's by D's and T's by A's. If we operate like that on the terms between parenthesis of i M 
in Eq. (2.9), we obtain after a tedious but quite straightforward calculation the following 
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result, has a term in one in D p £ a and one that contains the derivatives of the Lie 
derivatives of the background metric or of *** 



z pcj = \£ i g pa =D {p i rj) . (2.11) 

Here £ CT = <7 CTM £ M - Indices will always been displaced with the background metric , 
never with g pv . Thus has this form 

i M = (l? + ^9 pa R P *W + if) T + °» pa V P to + e M - (2.12) 
The undefined symbols in Eq. (2.12) satisfy the following equalities, 



2ni p = Q pa \{ A A x A' 1 + A M A A - 2 A*\ A A ^ - 8» ( A 71 

-I- n^ X ( A a A p — A a A p \ 



A\ - A\A A 

pa^rjX pX no 



(2.13) 



2«^ CT = (g^g™ + 9^ 9 pv ~ 9 pu 9 pa )Kx ~ (9 up 9 aX + 9™ 9 pX ~ 9 uX 9 pa )Kx ( 2 -14) 

and 

2/ce" = g pX d x z + g pa (d P z P(T - 2D p z£) , z = g p °z pcj = D x £\ (2-15) 

Had we applied Eq. (2.8) to C instead of C , we would have written everywhere g v 
instead of g^, from Eq. (2.8) up to Eq. (2.15). We would have found barred, conserved 
vector densities instead of that are as follows: 

¥ = (tit + r + (2.i6) 



The simpler form of Eq. (2.16) compared with Eq. (2.12) comes from the fact that A pa = 
and thus t„ = a^ pa = 0. Conserved vector densities for Co = C — C are thus obtained by 
subtracting Eq. (2.16) from Eq. (2.12). We find in this way, a conserved vector density 
relative to the background I M : 

I» = i p = ¥X + a^ a D p ^ + fj^, dj» = 0, (2.17) 

with 

¥ v = f» - f>? + ±-lP°R pa d p + t> (2.18) 



*** The presence of z pa comes from replacing second derivatives using the following iden- 



tity D pa e = W apv i v + 2D (p z^ - D P z pa . 
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and 



= ^ _ g M = J_ \l^d x z + l pa (d^z P(J - 2D p z$)\ . (2.19) 



For Killing vectors £ M , ff = and in Eq. (2.17), a^ pcr is equally zero. The remaining 
part of the cr-term contains the antisymmetric part a^ pa ^ to which we refereed in the 
beginning of this section. It plays the role of a (relative) helicity in linearized quantum 
gravity [30] and, see subsection 3(iii) below, is similar to the helicity in electromagnetism 
[31]. 

(iii) The superpotential. 

Had we not replaced (1/k)G p by T£ in Eq. (2.9) the conservation law <9 M i^ = would 
have remained an identity instead of becoming an equation. Thus i M must be equal to the 
divergence of an antisymmetric tensor density or "superpotential" [32] for any g^, g 
and i.e., there exists a i^ v = —i UfI such that i M = d^i^ . A superpotential is easily 
found by replacing in Eq. (2.9) £(T p l a by its expression in terms of D v £^g pcr : 

£^ pa = X -g^ {D p £ i9va + D a £ i9vp - D v £^g P ,) • (2.20) 

The conserved current itself 1^ = i^ — i^ is thus also equal to the divergence of a superpo- 
tential I pv = i^ u — i^ u . This superpotential is [see [6], on flat backgrounds see [16]; notice 
that flatness makes no difference in Eq. (2.21)] 

I*" = I - D^i"!) + -i [p k u] = + -i^k u \ I p = dj^. (2.21) 

K K K 

In Eq. (2.21) K^ u may be called the "relative Komar superpotential", relative to the 
background because (1/ k)D^^ u \ obtained with the Hilbert Lagrangian, is known as (half) 
the Komar [13] superpotential. I^ v is linear in and its first derivatives and, using Eqs. 
(2.4) and (2.10), can be written as follows 

j»u = IfA[^ Ar ] + p^^x with p^ x = j-g Xp D a (g p ^g^) . (2.22) 

where g^ is the inverse of g pv . The tensor density x is Freud's [33] superpotential on 
a curved background which has been written in this form already by Cornish [34] (more 
on this in section 5). The great advantage of a superpotential is, as we said, to replace 
volume integrals of complicated vectors by surface integrals of relatively simple tensors. 

(iv) How does I p depend on a divergence in the Lagrangian density? 

We found so far that Ncether's method applied to the perturbed Lagrangian density, 
Eq. (1.1) with Eq. (1.2), generates with every smooth vector £ M a conserved vector density 
fr, Eq. (2.17), the divergence of a superpotential fr", Eq. (2.21) or Eq. (2.22), and both 
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depend on the divergence <9 M /c M . The contribution of the divergence to our superpotential 
is apparent in (2.21). What is the contribution of the divergence to the different parts 
of the conserved vector density: the energy-momentum tensor density #^ and the helicity 
tensor density 6" MP<J ? To find this out let us write the divergence of the second term of Eq. 
(2.21) in a form similar to Eq. (2.17) with a term in £ M and one in D p ^ a and valid for any 
fc": 

d v ( = \d x (o^ a] ) T + -ti^D^. (2.23) 

The factors of and D p ^ a represent the respective contributions in to and a^ pcr . 
There is no contribution to If, in accordance with Eq. (2.1), we indicate by a prime 
the parts of those tensors that are independent of the divergence we may write, see Eqs. 
(2.17) and (2.21), 

K = O'l 1 + Y K (^Dxk x - D v k») , (2.24a) 

& H P a = & >PP* + fy-gUa _ %»gPA . (2.246) 
ZK V / 

A summary of the results obtained in this section together with comments has already 
been given in subsection (i). 
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3. Belinfante's method and Rosenfeld's identities 

(i) Motivations and summary of results. 

In this section we modify the conserved vector densities J M obtained in Eq. (2.17) 
using Belinfante's [20] trick. In classical field theory on a Minkowski spacetime there is a 
problem with the canonical energy-momentum tensor, the equivalent of = Q p L g pv here. 
The tensor is divergenceless but not symmetrical and therefore in Minkowski coordinates 
the angular momentum tensor is not divergenceless; the total angular momentum is not 
conserved. The reason is that it does not take account of the spin of the field. It is this 
situation that Belinfante did remedy by changing the canonical energy-momentum tensor 
in such a way that the total energy momentum would remain unchanged and even the 
local density would still remain the same in the appropriate gauge. Rosenfeld [27] found 
the same correction independently and in arbitrary coordinates. We shall use Rosenfeld's 
method in the next section for a different purpose. The Belinfante correction has been 
applied to gravity in general relativity on a flat background in Minkowski coordinates by 
Papapetrou [21]. He could then calculate the total angular momentum at spatial infinity 
and give a physical meaning to some of the irreducible coefficients in asymptotic solutions 
of Einstein's equations. 

Here we first apply the method to the electromagnetic field on a curved background. 
The quantities involved are familiar in electrodynamics and the results illustrate well the 
effect of Belinfante's modification. We then apply the same method to the conserved vector 
1^ of a perturbed gravitational field on an arbitrary background. We obtain in this way a 
new conserved vector density Z M , see Eq. (3.9), which generates a new energy tensor density 
Tjj 1 , Eq. (3.10) with Eq. (3.8). There is also a new vector density the analogue, in Eq. 
(2.17), of the sum a^P^Zpa + fj^ which is zero for Killing displacements; it is defined in Eq. 
(3.10) and written explicitly in Eq. (3.25). We construct also a new superpotential X^ v 
which is of great simplicity, see Eqs. (3.18) and (3.19). Parts of it are quite familiar in the 
weak field approximation. The new Z Ml/ discussed in some detail below in (iv) is remarkably 
linear in . Even more remarkable is that the new conserved vector is independent of any 
divergence that is added to the Hilbert Lagrangian. The energy-momentum tensor given in 
Eq. (3.33) with Eq. (3.34) is particularly interesting and its properties are brought forward 
by Rosenfeld's identities (3.27) and (3.28). Identity (3.28) together with Eq. (3.22) shows 
clearly that the energy-momentum tensor density is symmetrical only on Einstein 
space backgrounds defined by Eq. (3.31). Identity (3.27) shows that on such backgrounds 
the energy-momentum tensor density is also divergenceless. 

(ii) Belinfante's correction in electro-magnetism. 
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The familiar example of an electromagnetic field in empty space will help to see what 
Belinfante's addition does to the canonical energy tensor. Let y/^-gtf be the Lagrangian 
density on a flat background for simplicity. Thus R vpa = and in arbitrary coordinates: 

= -^V^r p 9" a F^F p(7 , = d„A v - d v A^. (3.1) 

We now repeat on y^—gC^ the operations performed on C from Eq. (2.7a) till Eq. (2.19). 
However, to avoid complications inherent to arbitrary £ M 's, we use only Killing vectors 
of the background for which z pa = but we stick to arbitrary coordinates. Then, in 
(daggered) notations similar to those of last section, we find that if Maxwell's equations 
hold D u F fJ,u = 0, the conserved vector is of the same form as 7^ in Eq. (2.17) without a 
yf". 

I^ = 0^ + a^d [p ^ h (3.2) 

Here o^ p,pa is antisymmetric in pa. The canonical energy-momentum tensor 6>V and the 
helicity tensor are respectively given by 

flty = --^ (f^ p D u A p - ^F prT F pa 5^j , a^ prT = -^-F^ p A a \ (3.3) 

The Belinfante modification consists in changing /1> to 

jl> = /|> + D V (S^ V %) (3.4) 

in which 

The modified current is now of the form 

j|> = rt^^, = -L (f^F p v + ^g pu F pa F pc ^j . (3.6) 

This is the familiar symmetrical, divergenceless, electro-magnetic field energy-momentum 
tensor. 

(iii) Belinfante's correction for the conserved vectors in general relativity. 

We now define, by analogy with Eq. (3.4), a new conserved vector density J M by 
adding to 7 M , see Eq. (2.17), a divergence of an anti-symmetric tensor density constructed 
with the anti-symmetric part a^lpa] Q f frvpa obtained from Eq. (2.14): 

X p = /> + d v (s^ p i P ^ = d v (i"" + S pup ^ = d v Z pv (3.7) 



47T 



(3.5) 
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with 

The divergence added to I M is the Belinfante addition in arbitrary coordinates. The 
vector density 1^ is linear in in z pa , see Eq. (2.11), and its derivatives D\z pa ; 1^ has 
no term in d[ p £ a ] anymore. The new conserved current is thus of the form 

j M = ff£* + z» = Q v ±w (3.9) 

with 

fi? = et + D p s» p „ = {d» pu + s» pu ) z p<J + if (3.io) 

while 

jpu = p v + gi*"P^ ( 3-11 ) 

It can be seen that if ^ is a Killing vector of the background £ M for which z pa = ?) M = 0, 
the conserved vector is simply given by 

j„ = f^ v = q v jpv, (3J 2 ) 

where = J^(£ = f) and = T^(£ = f). 

Now consider for a moment the tensor densities a^ pcr ^ which makes up S^ 9 . The 
contribution from any any to the cr-tensors have been singled out in Eq. (2.246) from 
which follows, in particular, that 

a v[p°\ = & >»[p<t] _ 2_gi4pif\. (3.13) 
2^c 

With this result we can rewrite S 111,(1 as a sum of a /c-independent part S ,p,up and a k 
contribution and we find that 

gP-v £ P _ cj/p,v _ Cjtfivp _ frip\pv\ _|_ frlp\pv\ _ fr'v[pp]_ (3.14) 

Notice that I p,v in Eq. (2.21) contains exactly the same £/c-term as S p,l p ^ p but with the 
opposite sign. Thus the new superpotential (3.11) is unchanged by <9 M /c M or any other 
divergence added to Cq for that matter. It depends only on the Hilbert Lagrangian for 
the perturbations R — R and the method used to generate it. If we take Eqs. (2.21) and 
(3.14) into account, we may, instead of Eq. (3.11), write 

J"" = K» v + S'» vp i p . (3.15) 

The divergence of this superpotential is equal to the new conserved vector J M . Conse- 
quently Tf and are independent of any divergence added to the Lagrangian, contrary 
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to the canonical energy tensor obtained by Ncether's method alone. A direct calculation, 
without using superpotentials gives of course the same results. Calculations are only more 
cumbersome. 

We shall now work out explicitly the formulas for and for . The explicit form 
of is obtained below with Rosenfeld's identities, 
(iv) The explicit form of the superpotential X^ v . 

The tensor a'^ pa ^ is obtained from Eq. (3.13) in terms of A's using Eq. (2.14) for 
and Eq. (2.4) for k^. The sum of the a"s that make up kS'^ V(T in Eq. (3.14) is 
found to give 

K ^ua = (^p^ _ A° pX g^) g u]x + (g ap g^ - g ap g x[A j + g pX ^ x g u]a - (3.16) 

On the other hand kK^ v defined in Eq. (2.21) can be written with the help of Eq. (2.10) 
in this form 

K K^ = lP^D p ^ + g p ^A^C- (3-17) 

The sum of Eq. (3.16) times £ CT and Eq. (3.17) gives the new superpotential in terms of 
A's. However, X tiV has a much nicer form in terms of D derivatives of the metric g Mi , or 
rather in terms of With Eq. (2.3) we can write X^ v in the following form 

= iiPl^D p C ] + V^ x £ x (3.18a) 
where the P-tensor plays now the role of the F tensor in (2.22): 

fipvp = _ V ^P = ]_D (j (gP»l»° - gP»l»° - g°»i»P + g™pp) . (3.186) 

2k V / 
Another telling and useful form of the superpotential is 

J/*" = - (j p ^D p ^ + ^DJ^ a - D^C) ■ (3-19) 

The superpotential (3.19) linear in and its derivatives is (not surprisingly) reminis- 
cent of various familiar expressions used in the literature and which are based on linear 
approximations or expansions to higher orders: 
(a) On a flat background: 

In Minkowski coordinates g^ u = n^ = diag(l, — 1, — 1, — 1). There are ten Killing 
vectors in the background. The four components of the Killing vectors of translations can 
be taken equal to ^ = 5^ with a = (0, 1, 2, 3) and D v ^ = 0; thus the four corresponding 
superpotentials are 

fii*™ = J_0 ct ^g™ + rfjji"* - r) a » g ^ _ rf*g va ) . (3.20) 
2k 
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This anti-symmetric tensor density p^ ua is sometimes quoted as the Papapetrou super- 
potential [21]. Considered as a linearized approximation p^ ua is the same quantity as 
Weinberg's [35] QP ua and Misner, Thorne and Wheeler's [36] dpH^ 13 . Equation (3.20) 
is identical with the linearized approximations of Freud's superpotential, Eq. (2.22), or 
of the Landau and Lifshitz superpotential written in arbitrary coordinates like in Cornish 
[34]. 

With the three spatial components for the Killing vectors of rotation which in Xp 
coordinates are given by ^ = (5^7]^ — 8pT] a y)X' r , becomes the Papapetrou superpo- 
tential for angular momentum as given in his original paper. 

(b) On Einstein space backgrounds: 

The linear approximation of Eq. (3.19) in arbitrary coordinates with Killing vectors £^ 
is equal to the Abbott and Deser [37] superpotential worked out on an Einstein space 
background [Eq. (3.31) below]. To obtain their full non-linear expression replace by 
-^gH^ u where H*" = g^g va H pa with H pa defined in Eq. (A. 13). 

(c) It is important to note that the correction to namely S^ up ^ p is homogeneous of 
order two in \p v and its derivatives. This is easily seen with Eq. (2.14) because 

2w"W = (l^Pg^ - g^PF^ A^ A - (l^Pg^ x - g^F^ A^ A . (3.21) 

Thus IP" U and X^ v are equal in the linear approximation as well on arbitrary backgrounds. 
From this follows [see KBL97] that Tp u provides the correct energy and linear momentum 
at spatial infinity. It gives also correctly the Bondi [15] - Sachs [38] energy and linear 
momentum at null infinity (see [39]). Details and proofs are given in appendix where we 
show also that the Abbott and Deser superpotential does not give the Bondi- Sachs linear 
momentum. 

(v) The new Z,p vector density. 

This vector density plays a role in conservation laws that are not associated with 
Killing vectors as in the examples mentioned in the introduction and in the cosmological 
applications of section 4. The vector density has also an important task in Rosenfeld's 
identities (see below). Z,p is defined in Eq. (3.10). It must be noted in Eq. (3.9) that the 
factor of D p ^ a is indeed symmetrical in pa: This follows from the definition (3.8) of SP ap . 
We need a new symbol for that factor which plays a role later on; let us set 

^gP-P<? _ ^gt^P — frP-vp _j_ CjVP — ^(MP)^ _j_ fr(P-v)P _ fr(P<r)V — ftl(vp)o _j_ £'(H<7)P _ 0-'(pcr)M_ 

(3.22) 

Thus 

2,p = *SP pa z pa + rf (3.23) 
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with f}^ defined by Eq. (2.19). An interesting form of *S ppa comes from using Eq. (2.14) 
for «j(mpV : 

*S^ pa = —D u (2l^Pg a > - g^l pa - l pv g pa ^ . (3.24) 
Plugging Eqs. (3.24) and (2.19) into Z p gives the following form "anti-symmetric" in z <-> / 

2kZ p = 2 [z pa D p l p - i pa D p z^ - (z p(J D*l p ° - l pa D^z p ^ + (l» v D v z - zDJ^ . (3.25) 

(v) The Rosenfeld identities. 

The conservation law <9 M T p = which holds for any smooth vector ^ contains deriva- 
tives of £ M of an order as high as 3. Thus with the help of Eq. (3.9), d p l p = can be 
written in the form 

d p l p = d p (f p n + d p Z p = Ke + ^D^ + ^ p \D {pa) ^ + ^%D {wa) C = 0. (3.26) 

This identity holds for arbitrary smooth £'s. Therefore all the properly symmetrized /3's 
must be identically zero, /?„ = = (3^ pa J = (3^ P l) = 0. This is the way in which 
Rosenfeld [27] obtained a set of identities and found the Belinfante correction. TJ 1 appears 
obviously in (3^ = and D^T^ in /3 U = 0. The remaining identities are independent of 
the energy tensor. In Rosenfeld's work, which is in arbitrary coordinates but on a flat 
background, f3 p v = was the most interesting relation. It related the canonical energy 
momentum generated by Ncether's method to the symmetric and divergenceless energy- 
momentum needed in Einstein's equations. The divergencelessness was shown by (3 V = 0. 

Here, however, the background is not flat and Rosenfeld's identities give different 
and interesting results. It is obvious that the calculations of the /?'s asks for a lot of 
rearrangements of factors that come exclusively from d p Z^. This is somewhat tedious but 
really straightforward. The resulting identities have in the end a nice form: 

A, = Dfi* - ^D^R pa = 0, (3.27) 

fo v = p p p gp» = f + D p (*S p P u ) - -l pp R u p = 0, (3.28) 
= *S ipa) " + D^ ppav = 0, (3.29) 

^P°> = with fro*!* = fr°pv = (l»iPg°> _ IPOgVA . (3.30) 

If we remember that *S ppu is symmetrical in /iv, see Eq. (3.24), we see from Eq. (3.28) 
that f pv = f vp if and only if 

=0 - R^ = -kg pv (3.31) 
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where A is necessarily a constant. Thus the energy-momentum tensor is symmetrical on 
backgrounds that are Einstein spaces in the sense of A.Z. Petrov [24]. de Sitter and Ein- 
stein's cosmological spacetimes belong to this class but in general Friedmann-Robertson- 
Walker spacetimes do not. This is why our formalism with arbitrary £ M 's is precisely good 
in relativistic cosmology and why may be important in that case. 

When the background is an Einstein space, Eq. (3.27) shows that is also diver- 
genceless. Identity (3.28) provides an interesting form for T 7 ^; using Eq. (3.24) in Eq. 
(3.28) we can write f» u like this: 

fiMu = _D p fip^ + —ip[»B»\ + J_[/p(m££) _ \p°W p v\. (3.32) 
2k 2k 

On a flat background we see, looking at Eq. (3.20), that T Ml/ is the second order deriva- 
tive of a tensor which was also obtained by Papapetrou [21]. Thus Eq. (3.32) is the 
generalization of his equation to curved backgrounds. 

The new energy-momentum tensor can be calculated from Eq. (3.10) with the use 
of Eqs. (2.13), (2.14), (2.18) and (3.8). It contains three types of terms, very similar to 
those of 9^ u : a symmetric matter energy momentum of the perturbations, a symmetric 
field energy-momentum tensor f^ u = t v ^ and two non-derivative couplings to the Ricci 
tensor of the background, the second of which being antisymmetrical: 

fi*» = {f^g^P - fw) + + —i pa R pa p iV + -l x ^R^. (3.33) 

2k k 

The field energy-momentum tensor density is the following terrifying homogeneous quadra- 
tic form in Z^", their first and second order derivatives: 

Kf ^ = i ^g P v _ r »fp^ D^ x pX + (j p °g x ^ - gP a l X ^) D a 

+ gT (l^AjxA^ + $ A "A<jA£ + A^A^^ - 2Aj A A<£$">") (3.34) 



+ 9 Xv 



\q^A a ,A p + (A* A { ^ - A1 A (At - A? A^) q v)p 

2i> pX^arj ' I per Xrj Xrj " per '-^Xp'-^Tja J i) 



On Ricci flat backgrounds T pu , Eq. (3.33), reduces to the expression found by Grishchuk 
et al [40]. 

(vii) Second derivatives in the energy tensor? 

Second derivatives of g^ v appear in the field energy tensor. This needs some comments. 
The canonical field energy see Eq. (2.13), is quadratic in first order derivatives. Thus 
that tensor density depends certainly on initial conditions only. This is the normal behavior 
of a conserved quantity. 
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Notice that the volume integral of the new conserved vector density is equal to a 
surface integral of X^ v in which there are no more than first order derivatives. This is also 
a suitable result. 

Consider now the local quantities X M themselves. Suppose that initial conditions are 
defined on a hypersurface at a given time coordinate x° = 0. Initial conditions are the 
metric components and their time derivatives (modulo Einstein's constraints and gauge 
freedom). The coordinate density is equal to J° = 1° + dk(S oka with k = 1, 2, 3. Only 
spatial derivatives of S appear in X° because S is anti-symmetric in the first two indices. 
Thus since 1° and 5 ,0fc<J ^ (J contain only first order time derivatives, X° itself contains only 
first order time derivatives and therefore even the energy-momentum density contains no 
more than first order time derivatives. The Belinfante correction adds only second order 
spatial derivatives. It is sometimes required, like in KBL97 that the energy-momentum 
tensor should contain no more than first order derivatives. This requirement may be 
unnecessarily restrictive. 

(viii) Conservation laws directly from Einstein's equations ? 

There is no difficulty to rebuild Einstein's equations G% = kT p from the conservation 
laws d v Z^ v = J M . To do this we rewrite given by (3.19) as follows 

K ±^ = -D^r^P + + lP^D p ^ (3.35) 

where Q v = D p l pu . It is useful to keep track of Q v because Q v = is the familiar 
(generalization of the) well known De Donder gauge condition. Thus, if we take the 
divergence of Eq. (3.35) and remember that d v T pu = X p we find, by replacing X M with 
T p ^ u + Z^ and by its expression given in Eq. (3.33) that 

2nd v X^ = (D p D p iP u - 2D { p<3 u) + DpQPgV + 2l p ^R u f } - 2l pa R^ p ^ & + 2kZ» 



2kJP = 



2k (f^g^P - f* v + f + 2\p^R v ) + l pa R p(T gP v 



i v +2nZP. 

(3.36) 

We can remove 2KZ 11 and 2i p ^R V p^ v from both sides of Eq. (3.36). The remaining ho- 
mogeneous linear expression in £„ is true for any £„. The factors of £„ on both sides of 
the equality must thus be equal. We are left with a set of equations that are of course 
Einstein's equations in which the left hand side contains all terms that are linear in 

D p D p iP u - 2D { pg u) + DpQPg^ - l pa R pa gP v - 2l pCT Rp pa u = 2k (f^g v)p - + f ^) . 

(3.37) 

On a flat background, in Minkowski coordinates these are Einstein's equations as they 
were written down by Papapetrou [21]. Equations (3.37) have also been given in this form 
by Grishchuk et al [40] on a Ricci flat background R pa = 0. 
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The linearized approximation on a flat background with the De Donder gauge condi- 
tion is readily recognized as the gravitational wave equations written in arbitrary coordi- 
nates: 

D p D p (^g^) = 2nf^. (3.38) 

Equation (3.37) is an interesting form of Einstein's equations from which we could have 
constructed our conservation laws d v X^ v = J M . But who has thought that by adding 
1kZ» + 2lP^Rp ] £ u on both sides of Eq. (3.37) we would get d v T* v = ±n 

4. Conformal Killing vectors, conservation laws and integral constraints 
in cosmology 

(i) Motivations and summary of results. 

Here we illustrate the theory developed in the previous sections with some applications 
in theoretical cosmology. Conservation laws have been used previously in cosmology (see 
introduction) and the following examples give potentially useful new formulas. 

We start by considering Friedmann-Robertson- Walker backgrounds with their 15 con- 
formal Killing vectors. We take the metric of the backgrounds in the form given by Eq. 
(4.1) and with Eq. (4.2) in terms of a conformal time rj. In these coordinates the conformal 
Killing vectors satisfy the 10 equations given in Eq. (4.6). A set of 15 linearly indepen- 
dent solutions of those equations is given in Eqs. (4.7) and (4.8) for spacetimes with flat 
spacelike sections r\ = (k = 0) and in Eqs. (4.7) and (4.10) for curved spacelike sections 
(k = ±1). Conserved quantities and integral constraints at a given time r] for perturba- 
tions of the background will be obtained from the time components of the conserved vector 
densities X° and the corresponding superpotential components X 01 (/ = 1, 2, 3). These com- 
ponents can be calculated for the 15 conformal Killing vectors and the general formulas 
for that are given by Eqs. (4.18a) or (4.186) and (4.19). Equation (4.18a) provide an 
expression for perturbations that may even be large. The rest of this section deals mainly 
with applications of these formulas. 

The 15 integrands *X and *X l defined in Eq. (4.20) are calculated explicitly. Notice 
that Eqs. (4.18) and (4.19) contain components of the conformal Killing vectors and 
their spatial derivatives. They do not contain the time derivatives of these components. 
We may thus use *X's and *X h s obtained from linear combinations of conformal Killing 
vectors with time dependent coefficients. The reason for doing this is that we find in this 
way 15 relatively simple integrands. They are written in Eqs. (4.24) for k = and in 
Eqs. (4.24) for k = ±1. Notice that linear combinations of conformal Killing vectors are 
not conformal Killing vectors anymore. Nevertheless with those non conformal Killing 
vectors we have obtained interesting "conservation laws". We show, for instance, that 
Traschen's integral constraint vectors [1] are equal to time dependent linear combinations 
of conformal Killing vectors. Formulae are explicitly given in Eq. (4.25) for k = and Eq. 
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(4.26) for k = ±1. We also show that if we equate to zero the uniform Hubble expansion 
rate Q [26] defined in Eq. (4.17), 14 of the 15 integral constraints take the particularly 
simple form shown in Eq. (1.10). The lone exception is associated with time translations 
(k = ±1) or time accelerations (k = 0). These exceptions are precisely those conservation 
laws that interested Uzan et al [8] (see also [7]) in which an unexpected field contribution 
in the conservation law of "energy" is inevitable. We also give a non trivial illustration of a 
spacetime that is asymptotically Schwarzschild-de Sitter with k = 0. In this case we obtain 
13 global integral constraints with no boundary contributions and two non-zero constants 
of motion for perturbations that may be large near the source but weak at infinity. The 
results, which may be new, are given in Eq. (4.34). For perturbations that are small 
everywhere, the formulas are given by Eq. (4.35). 

(ii) Friedmann-Robertson- Walker spacetimes and their conformal Killing vectors. 

We write the background metric ds 2 in dimensionless coordinates = (x° = n, x k ) 
with k, I, m = 1, 2, 3 for which the symmetrical role of x k is apparent: 

ds 2 = g^dx^dx" = a 2 (dr] 2 - f k idx k dx l ) = a 2 e llv dx il dx v ', (4.1) 

a{rj) is the scale factor and f kl and / = det(fki) are respectively given by 

fu = Ski + k \ X y** n , f kl =5 kl -kx k x\ /= r ^ 5 , (4-2) 
k = or ±1 and r 2 = 5kix k x l . The non-zero Christoffel symbols of the metric (4.1) are 
r° = d, T° kl = df kl , fS = d(JI», t^ = kx m f kh (4.3) 

a is the dimensionless conformal Hubble "constant" 

1 da , AA . 
a=- — . 4.4 
a an 

In these notations the non zero components of the Einstein tensor are respectively 

G° = \(k + a 2 ) = kT° , GT = \(k + a 2 + 2d a)5r = KTT, (4.5) 

Ct (X 

The conformal Killing vectors are the 15 linearly independent solutions £ M = £ fc ) of 
Eq. (1.8). These equations are independent of the conformal factor and can be written in 
3-dimensional notations as follows: 

d G e = \v k t\ d e = v^°, v( k e } = f ki o G e- (4.6) 
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where Vfc is a 3-covariant derivative for the fki metric, V = / V/, and the first equation 
equals one third of the trace of the last one. We found the solutions of Eq. (4.6) as follows, 
(a) The group of conformal transformations in Minkowski coordinates is explicitly 
given in [25]; infinitesimal transformations provide the conformal Killing vectors of the 
flat spacetime in coordinates, (b) The metric e Mi , is conformal to the Minkowski 
metric 77^, i.e. in appropriate coordinates e Mi , = 2 ?? Mi , globally and the conformal 
Killing vector components are the same since they do not dependent on O; we had thus 
only to take the components of £ M in Minkowski coordinates and transform them into our 
coordinates x^\ this is easily calculated from the explicit global coordinate transformations 
given in [41]. (c) One can in the end verify that the results satisfy indeed Eq. (4.6). Here 
are the results. 

There are 7 conformal Killing vectors which can be written in compact form for 
every value of k; these are the conformal Killing vectors of time accelerations (t), space 
translations (s a , a = 1, 2, 3) and space rotations (r a ): 

= sZ = 6Zy/l-kr 2 , rZ = 8» k e kal x l . (4.7) 

The other 8 conformal Killing vectors are somewhat different for k = and for k = ±1. 

For k = 0, the Lorentz boosts (l a ), dilatation (d), time acceleration (ao) and space 
accelerations (a a ) (the last two have been studied in [25]) are respectively given by 

l a = (l° a = x\ l k a = V 5 k ), d=(d° = r,, d k = x fc ), 

a = (ajj = v 2 + r\ a k = 2 V x k ), a a = (a° a = 2 V x a , a k = 2x k x a + [rf - r 2 ]5 k ). (4.8) 

(k = 0) 

For k = ±1 the 8 vectors can be written in a more compact form in terms of the 
column matrix 

j3 = ( ) = ( ^ ] with a = sin?] (k = 1) or a = smht] (k = —1). (4.9) 
\p. J \o aj 

What in flat spacetime corresponds to dilatation and time acceleration can be written as 
a single combination (5); the same is true of what correspond to the 3 Lorentz boosts and 
3 space accelerations (A a ) 

(6° = Py/l - kr 2 , 5 k = d py/l ~ kr 2 x k ), 

(fc = ±l) (4.10) 

= (\° a = d (3x a , X k =(3f ka )- 
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Notice that for k = we can take a = rj and apply Eq. (4.10) to that case as well. Then 
P = (i) and 5 = j while A a = )• These 4 vectors are not the same as those given 
in Eq. (4.8); only two of them are in that group d and l a , the other two t and s a belong 
to the group of Eq. (4.7). 

15 conformal Killing vectors are given by any linear combination (that are linearly 
independent of course) with constant coefficients of Eqs. (4.7) and (4.8) for k = 0, Eqs. 
(4.7) and (4.10) for k = ±1. 

Conformal Killing vectors and their space derivatives appear in the zero component 
of the conserved vector densities X° [see Eq. (4.19) below] in two combinations \z and a 
new one that we denote by y: 



4^ 



l Dp ^ = ±V k e + ae, y^{k + a 2 -d a)e + \(d G z-az) = \V 2 e + ke. (4.11) 



Most of the y's are zero. The non-zero one's are y(t) = k and y(ao) = 2. 

There are 3-antisymmetric tensors V' fc £^ which appear in the I 01 components of the 
superpotential, see Eq. (4.18a) or Eq. (4.186). The tensors are not zero for the following 
three conformal Killing vectors, 

V^ k s l l = -2kx^s l l V^r 1 } = e akl - 2kxK 1 }, (k = 0, ±1) 
V [k a l l=A5^x l \ (k = 0) 

(iii) Superpotentials and conserved vectors for small perturbations. 

We denote the perturbed metric components g^ u by + h^ v . Some authors [42] 
prefer to use the "conformal perturbations" and write g^ u = a 2 (e^ u + h^ v ). Thus, 

ds 2 = (g^ + h^)dx^dx v = a 2 (e llu + ~h^)dx»dx u . (4.13) 

We shall mainly use h^ u ; h^ v seems to be preferable in 4-covariant perturbation calcula- 
tions. In a 1 + 3 splitting, the 10 components of the perturbations are /too? hoi, h^i- We 
shall not displace the 0-indices up or down. The kl indices will be displaced with the fj-i 
metric. Thus 

hT = r k h kh h mn = f mk r r h kh ~K = r i h ol . (4.i4) 

There are simple relations between and h pu : = ^oo, h® = hoi and hf = —hf. 
The tensor density that enters in the superpotential (3.18) is related to h^ v in the linear 
approximation as follows 

= v^(-r p r a + \r v 9 pa )h P a = v7(-e w e w + § e ^o v- (4.15) 

Notice that Eq. (3.18) in terms of is the same for large or small perturbations and can 
be used for calculations at higher order of approximations. In this section, we are interested 
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only in small perturbations. There are two 3-tensors that deserve special notations because 
they appear as building blocks of the components of I 01 : 



(4.16a) 

In the linear approximation where Eq. (4.15) holds the quantities in Eq. (4.16a) reduce to 



q? = k?- sr~K, qt = (2d/ioo - v n h n)sr + v m h i - d oq r- 



(4.166) 



We also define by a special symbol Q the perturbed trace of the external curvature of 
the hypersurface n = const which appear in the zero component of the conserved vectors. 
Thus, if n M is the unit normal vector to that hypersurface, 



Q = - {-D^) = fd/ioo + \doK - V„/# 



(4.17) 



Q = is the "uniform Hubble expansion" gauge condition which was introduced by 
Bardeen [26]. 

We have now all the elements needed to calculate the conserved vectors and superpo- 
tentials with small perturbations for the 15 conformal Killing vectors and to write them 
down in a compact form. We are particularly interested in integral constraints over volumes 
at a constant time with spherical boundaries. For this we need only the zero components 
of the conserved vector: X° = dil 01 . Let us write first X 01 which is define by Eq. (3.18). 
We obtain after painful but straightforward calculations the following expression for the 
superpotential components which are valid for large perturbations, using Eq. (4.16a) 



-9 



2na 2 



(2ai 01 - V k q l k )e + q l k V h e + Q^* + '"Jv^ 



(4.18a) 



For small perturbations Eq. (4.18a) reduces with the help of Eq. (4.15) to 



J' 



0/ 



'-9 



2ko 2 



(2dh l - V fc gi)£° + q l k V k e + + h ok ^ k e ] 



(4.186) 



The linearized expression for X° is much simpler than it appears in Eq. (3.33) because 
T nis _ q j n terms of the energy- momentum perturbations 5Tjf = T£ — rather than 
perturbations of densities we find that 



'-9 



K,a z 



(4.19) 
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Equation (4.19) suggests that it makes sense to transfer V n {\zhQ ) from the left to the right 
hand side in X° = d{I 01 and to rewrite it in the following "renormalized" form appropriate 
to a 3-dimensional formalism: 



2 2 



*J = V;(*T Z ), *J=-^=2°-V l (J^), *2* = ^=X W - Jffcj,. (4.20) 

The zero indices are no more appropriate because the stared quantities are not the compo- 
nents of conserved vectors or superpotentials anymore. With the definitions in Eq. (4.20) 
we can write instead of Eq. (4.19) 

*J = U o e + U k e - yhl + \zQ, n M = Ka 2 5Tl (4.21) 

and instead of (4.186) 

*z* = \ (-vv*e° + iv*e° + - |^v*e* + ^v^) . (4.22) 

Integrating *X = V; over a sphere (r = const) at constant time 77 = const we obtain 

/ *I^/fd 3 x = r 2 ^/f <fs k i (*l k ) —sin{6)d6d(j). (4.23) 

We now give the list of the 15 *X's and their associated *X l, s. Some linear combinations 
with n dependent factors have greater simplicity. Such combinations break of course the 
group character of the algebra of globally conserved quantities but here we are interested 
in integral constraints at a given time 77 for which the group properties of our currents 
are not important here. We shall keep trace however of the corresponding 77-dependent 
combinations of conformal Killing vectors and use special symbols for *T and *l l that 
reminds us of their origin. For instance *Z(t) is denoted by T, *l(s a ) by S a and so on... 
. Thus for k = and k = ±1 we have the following *J's and *l h s with *J = V; 

t->T = ViT', T = U + 2dQ-kh^ = T -kh^ T l = -\V k q{, (4.24/) 



s a - S a = VjSi, S a = lWl-fcr 2 , ^^igiVl-^-^ofeAS, (4.24//) 

r a ^^ a = V/<, K a = Il k e kan x n , = ±(Q l k e kan x n + ~h ok e akl ) - kh ok xK 1 }. 

(4.24 ///) 

The following 8 quantities in which appears T defined in Eq. (4.24 I) are for k = only: 
l a - V s a - C a = Vj4, C a = T x a , 4 = -i V fe ^ a + (4.24 /F) 

d - 7/t -> D = V^', P = 5 a x a + 2Q, V 1 = \Q l k x k -h°\ (4.24 y) 
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SLo + v h-2 v d^A = V l A , Ao = %r 2 -2hl, A l = -\V k q kl r 2 + q l k x k , (4.24 VI) 
a a + r] 2 s a - 2r]l a -> 

A a = ViA l a , A a = 2Vx a - S a r 2 , A l a = Q l k (x k x a - \r 2 5 k ) + 2h ok (<f ^ - 5 kl x a ) . 

(4.24 VII) 

The next 8 linear combinations of conformal Killing vectors are for k = ±1 only. In those 
formulas the expressions like A a (/3) and \ a (dof3) represent the factors of (3 and 8q(3 in the 
conservation law associated with A . Thus 

X a (d p) - l a = V,4, £ a = T x a , C l a = -\V k q l k x a + \q a , (4.24 /V) 

<5(do/3) -> P = VjP', P = S a x a +2Q\/l - kr 2 , V 1 = {\Q l k x k -h l Q ) y/l - kr 2 , (4.24 V) 

6(J3) ^A = ViA l , A = T y/l-kr 2 , A l = -\{V k q kl + kq l k x k )^l - kr 2 , (4.24 VI) 

X a (P)^A a = ViA l a , A a = U k f ka -2kx a Q, A l a = \Q l a r} + kh l x a . (4.24 VII) 

We notice that V(k = 0) = V and V l (k = 0) = V 1 . Also £ a (Jfe = 0) = C a and C l a (k = 
0) = However, A (k = 0) ^ A and = 0) ^ A l as well as A a (k = 0) ^ A a . and 

(ii) Analysis of these results. 

T, S a , 1Z a and £ a , X>, Aq, A a for /c = or £ a , £>, „4 7 for = ±1 contain three 
types of terms: linear combinations of a 2 n5T® = 11^, the perturbation of the "uniform 
Hubble expansion" Q and the trace of the perturbation of the spatial components of the 
metric h™. 

S a and IZ a are homogeneous in n M . Besides these 6 quantities there are 4 additional 
linear combinations with ry-dependent coefficients that are homogeneous in n M . For k = 

a- 1 L a -lA a = d- 1 U x a -U k (x k x a -^ ka r 2 )^U,V a ^ 

, " (k = 0) (4.25) 

a^r - 1? = d _1 n - n k x k = n M y M . 

For fc = ±1 

a- 1 £ a + kA a = d- 1 U x a + kU k f ka = U fl V a ^ 

_ (k = ±1) (4.26) 

a^Ao -V = (a _1 n - n fc a; fc ) VI - At 2 = n M Vtf\ 

The 7's and V^'s are Traschen's [1] "integral constraint vectors" . Her vectors are thus linear 
combinations of conformal Killing vectors with time dependent coefficients. In particular 
the V a 's are linear combinations of A a (<9o/3)'s and A a (/?)'s while V is a combination of 

6(doP)'8 and 6(J3). 



29 



It is however clear that if we take the uniform Hubble expansion gauge condition 

Q = 0, (4.27) 

then 14 of the 15 "conservation laws" have volume integrands that are linear and homoge- 
neous in I1 M . The exception is for the conformal time translation (k = ±1) or acceleration 
(fc = 0): 

T = U -k~hl, (fc = ±l); A = n r 2 -2^, (/c = 0). (4.28) 

Thus if Q = 0, the conformal Killing vectors provide 14 linearly independent expres- 
sions that are momenta of order 0, 1 and 2 and are given by surface integrals involving 
boundary values only. Such expressions can in principle be constructed directly from 
Einstein's constraint equations. The constructs are however far from obvious, 
(hi) Example: spacetimes that are asymptotically Schwarzschild-de Sitter (k = 0) . 

In this example, the background is a de Sitter spacetime with k = and perturbations 
far away from its sources appear to be spherically symmetrical. However, perturbations 
may be large at and near the sources. The asymptotic metric in our coordinates has been 
given by Bicak and Podolski [43] [their formula (48)]. Neglecting powers of m/r higher 
than one, their metric is as follows 

ds 2 = dT 2 -(e 2HT + 2F)dx 2 -(e 2HT -F) X 2 (de 2 + sin 2 (e)d<p 2 ), F = ^ X ~ 3 e - H - (4.29) 

and H = (A/3) 1 / 2 ; m and A are constants. We set 

V = ee- Hr , r = H X , e = ±1 (4.30) 
so that Eq. (4.29) takes no this form 

ds 2 = (Ht])- 2 ^ 2 + (-8 k i + h kl )dx k dx 1 }. (4.31) 
Thus, comparing to Eq. (4.13), we see that 
a=(H\ v \)-\ />oo = 0, ^=0, h k i = s\mH (s M - 35 ^x™5 ln x^ ^ 

Notice that Q = in this example and furthermore h™ = 0. To calculate *l l we need the 
q's and Q's defined in Eq. (4.166): q[ = h l k and Q l k = —(3/rj)h l k . With these elements we 
find that all 13 integrals that follow are zero: 

I Tyffd 3 x= I S a ^fd 3 x= I U a ^fd 3 x= I C a Jfd 3 x= I A aV r fd 3 x = 0. 

J oo J oo J oo J oo J oo 

(4.33) 
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The equalities constitute as many Traschen-like integral constraints. The 2 constants of 
motion that are not equal to zero are associated with dilatations and with time accelerations 
in Minkowski spacetime which is conformal to our background 




e 



L 



OO 



V^ffafx = HMc 2 , 




e 



L 



OO 



Aoy/fd 3 



X 



= iHMc 2 . 



(4.34) 



In these expressions M = mc 2 /G. Notice that the 15 integrands of Eqs. (4.33) and 
(4.34) do not have to be linearized. Near the sources perturbations may be large. If the 
perturbations are weak in the whole space, we may write instead of Eq. (4.34) 



These are integral constraints on 5T® and 8T®. 

5. Comments on the role of superpotentials in the theory of conservation laws 

(i) Motivations. 

Here we want to connect our work with past literature, give due credit to yet un- 
mentioned papers and make contact with some well known superpotentials or energy- 
momentum tensors that we have not yet encountered. 

(ii) On superpotentials in conservation laws today. 

Perhaps the single most important legacy of studies on conservation laws in general 
relativity is that conserved quantities in finite volumes can always be expressed as surface 
integrals on the boundary of the volume. Anti-symmetric tensor densities like X^ v dominate 
the scene to day in the literature, not X M . One great push in that direction was given by 
Penrose [14] who introduced the notion of "quasi-local" quantities which, in the weak field 
limit, reduce to ordinary conserved linear and angular momentum of the gravitational field 
in finite volumes. Many papers have been published on the subject in particular on quasi- 
local energy. Unfortunately, selection rules are few and no consensus exists. There is an 
interesting comparison of formulas in a paper by Berqvist [44] on the energy enclosed by 
the outer horizon of a Kerr black hole in which it is shown how six different formulas give 
five different results. There exists however a common point to those various definitions 
of quasi-local energy: they are not related by differential conservation laws to Einstein's 
equations [45]. In this instance, the present work points in a very different direction. 
We have not tried to make the connection with our own superpotential. The role and 
importance of superpotentials in field theory has been emphasized by Julia and Silva [28] 
who gave them an elegant and general mathematical basis. 

(iii) Connection with other superpotentials on a flat background. 




5T Q k x k ^]d 3 



x = HMc 2 , — / 5T°r 2 v r fd s x = \HMc 2 . 



T Joo 



(4.35) 
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Rosen [19] was the first to drive attention to the fact that the quadratic IT-Lagrangian, 
Eq. (2.6) used by Einstein to derive a conserved pseudo-tensor could be written in covariant 
form by introducing a second metric. This amount in practice to describe curved space- 
times with respect to a flat background. The mathematical basis of Rosen's approach is 
given in Lichnerowiscz [46]. It is thus no surprise that our formulation of conservation 
laws for perturbations of curved backgrounds connects nicely with well known conserva- 
tion laws in classical general relativity. This is what we want to show here. Let us go back 
for a moment to the "divergence dependent" conserved vector J M with Eq. (2.17), not 
X M . The Rosenfeld identities have been worked out in KBL97****. On a flat background, 
R vpa = and in arbitrary coordinates, the formulas of KBL97 are similar to Eqs. (3.27) 
- (3.30). In particular, 

= 0, 6» v = -D x a x ^, 6$ = T» + %, (5.1) 

t% is Einstein's energy- momentum tensor density as given by Rosen in arbitrary coordi- 
nates. On a flat background 0£ is thus the divergence of a tensor, not an anti-symmetric 
one and not a two index tensor but still one acting like a "superpotential" for volume 
integrals in Minkowski coordinates. a x ^ v is Tolman's [47] superpotential, apparantly the 
first of its kind in the literature. This superpotential is related to another famous one, 
Freud's [33] superpotential F XfJ,u is defined in arbitrary coordinates by Eq. (2.22): 

<r A ^ = F x ^ + —D p (g u H x ^) . (5.2) 

Since covariant derivatives on a flat spacetime are commutative, by taking the divergence 
of a x w and using its relation with 6^ v - see Eq. (5.1) - we obtain a similar relation 

= -D X F X ^. (5.3) 

The great simplicity of Freud's superpotential made it a successful quantity to calculate 
globally conserved quantities like the total energy at spatial infinity as well as at null 
infinity [48]. Our own energy tensor on a flat background satisfies similar relations. 

D lt tf = 0, f » v = -D x (*S XfM/ ) = -D X V X ^. (5.4) 

The difference between T^" and 9^ v following from Eq. (3.10) is exactly D p S p ^ v . It is 
interesting to notice the relation between the divergences of the *S and V tensors when 
backgrounds are not flat: 

d x v x ^ = D X (*s x n - y/^p* - y/ v ^ x - (5 - 5) 

**** In KBL97, the formulas are (2.50) to (2.52) but a " is missing on every symbol! 
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A most famous superpotential is that of Komar, (1/k)D^^. Its greatest quality is 
to be background independent. It is also useful in asymptotically flat spacetimes at spatial 
infinity. But it has some shortcomings which we already mentioned. There have been 
various corrections of that attractive covariant expression [49] which did not get rid of the 
anomalous factor 2 and had also other "defects" [45] . 

One intriguing superpotential is that of Arnowitt, Deser and Misner [50], especially 
for energy. It was original defined in a synchronous gauge, at least asymptotically, goo = 1, 
g k = 0. In that gauge, the surface integral at spatial infinity of the Komar tensor is 
zero. What remains then of the superpotential in Eq. (2.21) is the (Ok) -component of 
(1/ ' K)^k v ^ which reduce indeed to the ADM integrand at infinity as can easily be verified. 

The Landau and Lifshitz superpotential is as satisfactory as our own superpotential 
for calculating the total 4-linear momentum but it has the wrong weight and it is difficult to 
see how to connect it with the group of diffeomorphisms via Ncether's method on a curved 
background. The L-L complex has however been obtained recently from a variational 
principle by Babak and Grishchuk [51] on a flat background in arbitrary coordinates. 
Incidentally L-L's pseudo-tensor has one more drawback, not shared by Einstein's pseudo- 
tensor which was pointed out by Chandrasekhar and Ferrari [52]. Consider the weak 
field approximation of the total energy in a stationary spacetime. A variational principle 
applied to the total "Einstein Energy" leads to Einstein's linearized field equations. The 
"L-L Energy" provides incorrect equations. 
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Appendix on Globally Conserved Quantities 

(i) Object of this appendix. 

There exists quite a number of different conservation laws but few criteria to select 
among them. The main discriminating conditions are global conservation laws and the 
weak field limit which are to some extend related. The only non-ambiguously defined 
global quantities are the 4 components of the linear momentum P a at spatial and at null 
infinity for spacetimes with definite fall off conditions to asymptotic flatness. It is therefore 
important to show that our superpotential gives at least the correct total 4-momentum in 
those cases. In this task we shall be greatly helped by Eq. (3.11) which for Killing vector is 

= ji»> + s» v %. (A.l) 

It has been shown that the superpotential J^ v provides by itself the total 4-momentum P a 
both at spatial infinity [16] and at null infinity [53]. We must therefore show that S^ vp ^ p 
does not contribute to P a in both asymptotic directions. This is what we briefly indicate 
in this appendix. We shall show that in both asymptotic directions the "discrepancy" 

AP a = £ S^ {a)p dS^ = 0. (A2) 

S is the sphere at infinity and £( a )p 5 ( a = 0, 1, 2, 3) are the four Killing vectors of transla- 
tions in the asymptotically flat background. 

(ii) Stationary spacetimes and spatial infinity. 

Consider stationary solutions that fall off as follows in asymptotic Minkowski coordi- 
nates at spatial infinity x p = (x° = t, x k ): 

1 1 

9nu(x ) = r}^ + -u pv + 2 , dig^(x k ) = ^v iliv + 3 (A3) 

where r = \jY,(x k ) 2 and u pu , vi pu depend on the directions on the sphere at infinity. 
The background metric is r\ pv and the Killing vector of translations i( a ) u = r\ av . The 
discrepancy is therefore given by 

AP a = I S ola m r 2 sin6d6d(f>, n t = — . (A4) 

It takes some patience to make this calculation, using the asymptotic form of the metric 
(A. 3) in Eq. (3.8) with Eq. (3.21), but the calculation does not need clever tricks and we 
find indeed that AP a = for r — > oo. The metric does not have to fall off as fast as in Eq. 
(A. 3). AP a = under weaker fall off conditions has been studied in [54]. 
(in) Radiation superpotentials at null infinity. 
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For the Bondi-Sachs metric we use the Newman and Unti [55] representation with 
coordinates x x = (x° = it, x 1 = r, x 2 , x 3 ) in which the metric has the following form 



ds 2 = goodu 2 + 2dudr + 2g Lcludx L + gKhdx K dx L (A5) 
and the flat background has a 

ds 2 = du 2 + 2dudr - ^[(dx 2 ) 2 + (dx 3 )% P=\ + \[(x 2 ) 2 + (x 3 ) 2 ]. (A6) 

The asymptotic form of the metric depends on 5 independent real functions of (u, x 2 ,x 3 ), 
namely ip", ip' 2 , u' and a" . These notations are similar to those of Newman and Unti - 
without a zero index - who use complex functions; here the prime indicates the real part, 
two primes indicate the complex part of their complex functions. The asymptotic form of 
the metric in these notation is 

2ib' 

<7oo = l + — + 2 , (A7) 
r 

902 = -P 2 ld2(^ 2 )+d3(^)]+l^+0 2 , g 03 = P 2 [ds(^)-d 2 (^)] + 2 -^-+0 2 , (A8) 

// 2 / 2 / 

U T f(j f TU 

923 = -r-^ #22 = -^--^ + Oo, 933 = ~^ + ^ + O . (A9) 

with \a\ 2 = a' 2 + a" 2 . 

The Killing vector of translations in the flat background spacetime with the metric 
(A. 6) have the following components 

- - X m 

£(o) M = (1, 1, 0, 0), £(m) M = (0, -n m , -rd 2 n m , -rd 3 n m ) , n m = — — . (A10) 

With Eq. (A. 5) to Eq. (A. 10), we have the necessary elements to calculate Eq. (A. 2) at 
null infinity using Eq. (3.8) with Eq. (3.21). The calculation is even more tedious than 
before but still it is not difficult to show that indeed AP a = for r — > oo. 

The reader may remember that the loss of energy E per unit u-time is given by 

^ = -(87T)" 1 j[{d u a'f + (d u a") 2 ]p- 2 dx 2 dx 3 < 0. (All) 

This formula is one of the outstanding results of Bondi. 

The loss of energy obtained from the Abbot and Deser [37] superpotential is differ- 
ent and is not negative definite. Following the prescription indicated in section 3 this 
superpotential satisfy the following expression 

kJ^ d = -H p ^D p ^ + l P D {il R v } - l^D p H^ p (A 12) 
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where 



(A.13) 



If we use this superpotential to calculate the energy Ead and the corresponding energy 
loss, we find that 



This is not negative definite. The Abbott and Deser superpotential should be rejected on 
this basis. 
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